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GHAVSAT/GEOPAUSE COVARIANCE ANALYSIS 
INCLUDING GhOrOTENTIAL ALIASING 

David VV. Koch 


ABSTRACT 

NASA's ongoing Earth and Ocean Dynamics Applications Program 'EODAP) 
uses satellite technology and instrumentation to sense and monitor a wide range 
of natural phenomena. Among the proposed missions within this program is the 
GHAVSAT/GEOPAUSE which employs satellite-to-satellite tracking technology 
in an attempt to obtain improvements in the accuracies of the geopotential field 
and the global geoid. This report describes a conventional covariance analysis 
of that mission in which the uncertainties of approximately 200 parameters, 
including the geopotential coefficients to degree and order 12, are estimated 
over three different tracking intervals. The estimated orbital uncertainties for 
both GHAVSAT and GEOPAUSE reach levels more accurate than presently avail- 
able. The adjusted measurement bias errors approacli th»> mission goal. Survey 
errors in the low centimeter range are achieved after ten days of tracking. The 
ability of the mission to obtain accuracies of geopotential terms to (12, 12) one to 
two orders of magnitude superior to present accuracy levels is clearly shown. 

A unique feature of this report is that the aliasing structure of this (12, 12) 
field is examined. It is shown that uncertainties for unadjusted terms to (12,12) 
still exert a degrading effect upon the adjusted error of an arbitrarily selected 
term of lower degree and order. Penally, the distribution of the aliasing from 
the unestimated uncertainty of a particular high degree and order geopotential 
term upon tne errors of all remaining adjusted terms is listed in detail. 


GRAVSAT/GEOPAUSE COVARIANCE ANALYSIS 
INCLUDING GEOPOTENTIAL ALIASING 


EODAP AND GRAVSAT/GEOPAUSE 

Tlie Earth and Ocean Dynamics Program (EODAP) formerly EOPAP (Ref. 1), 
conceived and managed by NASA, applies sophisticated satellite technology to 
the continuous global sensing and monitoring of a wide range of diverse natural 
phenomena. Earth orbiting satellites, utilizing their complex onboard sensing 
instruments, provide the unprecedented capability for automated global data 
gathering. The EODAP program is specifically designed to exploit the 
inherent advantages of satellites and their instruments to achieve a better 
understanding of earth and ocean dynamics. Continuous and precise observations 
of the earth's various dynamic phenomena are necessary to develop the complex 
mathematical models of earthquakes. Accurate observations of the ocean's 
dynamics lead to increased understanding of the surface conditions and general 
ocean circulation. Such data can be utilized immediately by the shipping in- 
dustry. It needs this data to plan ship routes which can take advantage of pre- 
vailing currents with attendant savings in both time and fuel. Indeed, it is 
difficult at present to foresee the full potential of this monitoring of natural 
phenomena by satellites. 

Because of the variety and complexity of natural phenomena, the EODAP pro- 
gram necessarily involves many satellite missions, e.g. the GEOS-C, MAGSAT, 
SEASAT, GRAVSAT, and GEOPAUSE. GEOS-C employs a satellite-to-ocean 
radar altimeter with a one meter precision to determine the short- wavelength 
undulations of the oeean geoid to an equivalent vertical accuracy. SEASAT, a 
successor to GEOS-C, uses this same altimeter with a pro|x>sed 10 cm accuracy 
level to sense the sea surface, currents, and general ocean circulation. Such 
data will be used in the study of geoid undulations, wind stresses on surface water, 
location and mapping of boundary currents, and tides. GHAVSA'l' (Bel. If) is a pro- 
posed low altitude (~:f00 km) surface force compensated satellite used with a com- 
panion satellite at the same height or with a high altitude satellite such as 
GEOPAUSE or |X)ssibly a tracking and data relay satellite (TDRS). GIUVVSAT 
is an instrument which is expected to provide the means for a determination 
of the geopotential field to a 2.5” resolution and a global geoid to a 10 cm ac- 
curacy. An improved detennination of the field will have an immediate Ijene- 
ficial impact. In satellite orbit determination, the geo|X)tential inaccuracies 
are a well established error source, especially over long orbital arcs. ;\s an 
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example, Ref. 3 describes the degrading effect of these inaccuracies upon a 
multi -revolution GEOS-C orbit. Satellite tracking and the subsequent data 
processing could be performed more efficiently and economically with a more 
accurate field, since in a given situation less data would have to be processed 
to determine an orbit. 

GEOPAUSE (Ref. 4) is a high altitude tracking relay spacecraft in a pre- 
cisely determined orbit. Data from this satellite is expected to determine the 
orbits of lower satellites such as SEASAT and GRAVSAT to a 10 cm accuracy 
and also to establish the position and dynamics of points on the earth's surface. 
The latter includes determination of polar motion, earth rotation, tectonic plate 
motion, and solid earth tides. Furthermore, the data may also be used to esti- 
mate both reference station locations and radio tracking instrument biases to 
the centimeter level. 

The basic GEOPAUSE satellite system concept has been investigated pre- 
viously by Siry (Ref. 5). Therein he describes the potential contributions and 
characteristics of such a system. In addition, the utility of the system in geo- 
dynamics and global surveys of natural phenomena is discussed. Another report 
by Koch and Argentiero (Ref. 6) documented an initial mission analysis study of 
the combined GRA VSAT/GEOPAUSE concept oriented toward low degree and 
order geopotential uncertainty recovery. That particular report demonstrated 
that data from GRAVSAT/GEOPAUSE could be used to determine uncertainties 
in the spherical harmonic coefficients of the geopotential field from (0, 0) through 
(8, 6) inclusive to an accuracy one to two orders of magnitude over current 
levels. This present report follows and expands upon that study. In the present 
report the uncertainties of a larger (12, 12) field are recovered, additional 
perturbative forces are considered, and parameters such as survey errors 
and instrument biases are estimated. An important matter, usually not investi- 
gated, is the aliasing structure of such a larger geopotential field. The aliasing 
effect of the uncertainties of unadjusted geopotential terms upon the error of an 
arbitrarily selected adjusted term is considered herein. The present report 
also stresses the difference, if any, between uncertainties of various parameters 
when the tracking interval is varied over three, six, and ten days. 
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THE GRAVSAT/GEOPAUSE MISSION 

TTie proposed GRAVSAT/GEOPAUSE mission involves satellite-to-satellite 
tracking technology combined with a unique orbital geometry in an attempt to 
obtain improvements in the accuracies of the geopotential field and global geoid. 
Table 1 lists, in a general context, the rationale, description, and the applica- 
bility of the results for such a mission. 

Table 1 

The GRAVSAT/GEOPAUSE Mission 


I. RATIONALE 

Ihe present uncertainties in the geopotential field are a fundamental error 
sour . >3 in earth and ocean physics. Meeting stringent EODAP objectives 
will require a new approach to estimation of these accuracies. 

II. DESCRIPl^^N 

Hi-lo satellite-to-satellite tracking data used to extract field 
coefficients from low satellite orbital perturbations. 

(1) GRAVSAT 

Low altitude surface force con.pensated satellite 

(2) GEOPAUSE 

High altitude tracking relay satellite. 

III. APPLICABILITY OF RESULTS 

More accurate geopotential field leads to improvements in: 

(1) Orbit determination and prediction 

(2) Satellite data processing 

(3) Accuracy of earth and marine geoid. > 


The overall geometry of the mission is depicted in Figure 1. 'I'lie mission 
involves a low ilRAVSAT and a high GEOPAUSE satellite in coplanar, polar, and 
circular orbits. Their common orbital plane is normal to both the ecliptic :md 
equatorial planes. The perpendicularity of the plane to the ecli|)tic minimizes 
long term luni-solar perturbations u|)on the plane, thus insuring long term orbital 
' stability. A jwlar orbit guarantees complete sampling of the earth's gravity 
field by GRAVSAT. Geojwtential coefficients from such a global sample will 
achieve orthogonality and will not be degraded by aliasing or |X)or statistical 
imlependcnce. The orbit of the high GEOPAUSE is minimally |>erturl)ed In- 
existing inaccuracies in the gravits field. Conse(|uently GE( )P.\l’SE can lunction as a 
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GEOPAUSE 

h - 3.6 E.R., CIRCULAR AND POLAR ORBIT 
P - 14 HRS. 

LOCATED AT GEOPAUSE 
ACCURATELY DETERMINED ORBIT 


> - RANGE SUM RATE 

MEASUREMENTS 


GRAVSAT 

h - 300 KM. 

P» 90MIN. 

ORBIT CIRCULAR AND COPLANAR 
WITH GEOPAUSE 
ORBIT HIGHLY PERTURBED BY 
DRAG & GRAVITY 

SURFACE FORCE COMPENSATION (SFCI 
SYSTEM REMOVES DRAG PERTURBATION 


FiKure 1. GHA VSAT/GEOPAUSE Geometry 
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stable high altitude relay tracking station. The GRAVSAT orbit at a 
height of only 300 km will be highly perturbed by external forces, principally 
atmospheric drag and gravitational uncertainties, and to a lesser extent by im- 
plementation of surface force compensation (SFC). Two difficulties immediately 
arise. First, the orbital perturbations due to external forces must be separated 
as to their origin. Second, the perturbations due to implementation of SFC must 
be properly treated so that any lack of knowledge of this self gravitational force 
will not degrade the remaining estimates. 

The first difficulty is overcome by the unique construction of the actual 
GHAVSAT spacecraft combined with an onboard SFC. The spacecraft will 
consist of two concentric spheres. The outer sphere experiences perturbations 
caused by both gravitational and non-gravitational forces whereas the inner 
experiences ''olely the former. This imbalance of forces causes the outer 
sphere to be tux>gressively displaced toward the inner as the GRAVSAT pro- 
ceeds in orbM The sensitive instrumentation of the SFC system will sense 
such a displacement. When this displacement exceeds a prescribed distance, 
the SFC system actuates jets on the exterior of the outer sphere until both 
spheres are concentric again. Meanwhile, during the entire displacement and 
repositioning, the inner sphere has followed an orbit perturbed solely by gravity. 
Ihus the SFC essentially has nullified the orbital perturbations due to non- 
gravitational forces and has constrained both spheres of the GRAVSAT to follow 
an orbit perturbed solely by gravity. It is these orbital perturbations, caused 
only by gravity, that the tracking system senses and from which gravity coeffi- 
cient estimates can be extracted. 

The secx)nd difficulty involves the proper statistical treatment of the force 
caused by the mutual attraction between the masses of the inner and outer 
spheres of the GRAVSAT. This force is assumed constant of unknown magnitude 
and acting solely in the along track direction op|X)site to the direction of the 
velocity vector. 'Hierefore, it can be treated as a standard differential correc- 
tion type bias and its uncertainty estimated. 

A tracking network of six ground stations was chosen to track GRAVSAT 
throuf^ GEOPAUSE. Three properly placed stations in each hemisphere pro- 
vide continuous coverage of GEOPAUSE by at least one station. Such coverage 
insures no permanent spatial gaps in the sampling of the gravity field by 
GRAVSAT. Each station measures range sum rate with a one minute integration 
time over a ten day data interval. This interval constitutes one c*omplete 
sampling of the gravity field by GRAVSAT. 


The nominal mission parameters follow: 


(1) ORBITS 

Epoch: July 1, 1979, 0“, UT 

Geocentric True Earth Equator and Equinox of Date 



GRAVSAT 

GEOPAUSE 

a(km)* 

6,678.133 

29,431.213 

e 

0.0 

0.0 

i(deg) 

90.0 

90.0 

n (deg) 

90.0 

90.0 

a (deg) 

180.0 

180.0 

M(deg) 

180.0 

180.0 

h(km) 

300.0 

23,053.190(3.6 E.R.) 

P(hrs) 

1.5 

14.0 


/ 


(2) TRACKING STATIONS 

Northern Hemisphere: Guam, Madrid, and Rosman 
Southern Hemisphere: Canberra, Johannesburg, and Santiago 

(3) MEASUREMENTS AND UNCERTAINTIES 

Range sum rate measurements every minute over ten days with random 
0.2 mm/sec noise and fixed mm/sec bias. 

(4) ESTIMATED PARAMETERS 


NAME DIMENSION 

GRAVSAT epoch state 6 

GEO PA USE epoch state G 

Range sum rate bias G 

Station survey components 15 

Residual force upon GRAVSAT 1 

Sine and cosine components of 1G5 


spherical harmonic expansion 
of geopotential field through (12, 12) 


*a through M are orbital elerrtents; a is the semi-major axis, e is the eccentricity, i is the inclination, 
il is the right ascension, cj is the argument of prerigee, and M is the mean anomaly, h and P are the 
orbital height and period, respectively. 
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RESULTS 


Under the previous assumptions, a rovariance analysis simulation was per- 
formed using the NAP parameter estimation/NAPCOV covariance analysis 
computer programs (Ref. 9-12). Instead of performing an actual least squares 
reduction of simulated or real data, the simulation derived only the covariance 
matrix of such a process. Furthermore, the simulation assumed that perturba- 
tions in the data are linearly related to the perturbations of the systematic and 
dynamic error sources under consideration. The uncertainties for all parameters 
were estimated simultaneously. Solving for all parameters together should 
theoretically minimize the uncertainties of their estimates. However, this ap- 
proach ignores the effect of the uncertainties of the unadjusted parameters upon 
the errors of the estimated ones (aliasing) and consequently may give optimistic 
results. 

Tables 2 through 6 list the a priori and a posteriori uncertainties after the 
processing of three, six, and ten days of tracking data for all the parameters 
estimated during the numerical simulation. For convenience, both the individual 
and RSS* uncertainties for groups of parameters are listed where applicable. The 
correlation matrix between the GRAVSAT and GEOPAUSE epoch states after ten 
days of tracking is also listed. The first two tables concentrate on the adjusted 
errors of the GRAVSAT and GEOPAUSE epoch states and their associated correla- 
tion structure. Table 2 indicates a substantial overall reduction in uncertainties 
of each component of both state vectors from their rather conservative a priori 
values. This reduction of uncertainties demonstrates the ability of the GEOPAUSE 
to act as a tracking relay spacecraft for use in recovering low satellite orbits. 

For the GRAVSAT orbit the largest uncertainties are predictably out-of-plane 
(X direction). Since their two orbits are coplanar, GEOPAUSE can sense only 
in-plane (y and z directions) GRAVSAT orbital perturbations and hence no out- 
of-plane information about the GRAVSAT orbit is obtained. However, with increased 
tracking data even these out-of- plane errors tend to decrease. The a posteriori 
RSS position uncertainties drop approximately linearly with time whereas the 
corresmnding velocity errors drop by a factor of approximately twenty percent 
between tracking intervals :md eventually reach the 2 mm/s level at 10 days. 

Ik)ln position and velocity adjusted RSS uncertainties are almost entirely com- 
posed of out-of-plane errors. The remaining estimated in-plane uncertainties 
are relatively insensitive to additional tracking. 

For the GEOPAl’SE orbit, the out-of-plane and one in-plane (z direction) 
a posteriori ;x)sition uncertainties behave in a similar m:mner. Roth decrease 
slowly and eventually reach a U vel of alx)ut 05 cm alter ten days of tracking. The 
adjusted error of the remaining position coni|X)nent remains I'ssentially constant, 
riu- adjusted RSS position uncertainty decrt‘ases uniformly from 2.11 nulers 
alter three days ol tracking to Ixdow the meter level alter ten days. ’I he adjusted 
RSS velocity uncertainties remain below 0.25 mni/ s tlirougluxit. Aeeoiding to 

* RSS IS The sciuari* root of the sum of the squares 
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3, 6, and 10 Days of T»'acking 


Epoch States 

A Priori 
Uncertainty 

A Posteriori 
Uncertainty 

3 Days 

G Days 

10 Days 

X (m) 

200 

mmm 

8.14 

4.78 

Y (m) 

200 

■M 

.02 

.02 

Z (m) 

200 

.35 

.20 

.14 

RSS (m) 

34 G 


8.14 

4.78 

onAVSAT 

200 


13.32 

2.3G 

t (mm/s) 

200 

.41 

.23 

.17 

Z (mm/s) 

200 

.07 

.04 

.03 

RSS (mm/s) 


72.25 

13.32 

2.37 

:• (m) 


1.43 

.89 

.G5 

' (HI) 


.07 

.04 

.03 

Z (m) 

150 

1.55 

.87 

.G2 

RSS (m) 

2G0 

2.11 

1.25 

.90 

GEOPAUSE - r. - - “ 7 “ 

X (mm/s) 

150 

.13 

.07 

.05 

Y (mm/s) 

l.'^O 

.19 

.11 

.08 

(mm/s) 

150 

.01 

.00 

.00 

RSS (mm/s) 

2G0 

.23 

.13 

j .09 


published literature the present state-of-the-ari orbital position uncertainty is 
approximately 10 meters. Hesults indicate that data from GEOPAUSE tracking 
GRAVSAT can be processed to achieve such an orbital accuracy within six days. 

Table 3 lists the symmetric correlation matrix between the GRAVSAT and 
GEOPAUSE epoch states after ten days of tracking. In general, the correlations 
reveal a good separability between estimates which is reflected in the a posteriori 
uncertainties of the previous tab!.. The correlation coefficients conveniently 
separate into two distinct classes. The first consists of correlations with absolute 
values greater than 0.94 and the second with values less than 0.43. The correlation 
coefficients of large absolute value are circled within the table. These high cor- 
relations portend possible convergence i roblems during any future processing of 
real data from the mission. TTie highly correlated parameters are the /. and Y 
components from either state and GEOPAUSE Y with Indeed, two of these 
pairs, GEOPAUSE Y with Z and GEOPAUSE Z with Y have correlations with 
absolute value near unity. All remaining pairs of parameters are relatively 
weakly cxirrelated. 
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Table 3 

Correlation Matrix of GHAVSAT and GEOPAUSE States After 10 Days of Trackinj^ 


GRAVSAT 


X 

Y 

z 

X 

Y 

X 1.00 

.03 

-.13 

.42 

.13 


1.00 

-.02 

.09 

.03 

7 


1.00 

-.18 


<x 



1.00 

.17 





1.00 


Z 

X 


Sv SYMMETRIC 

w Y 
Z 


GEOPAUSE 

Z X Y Z X Y Z 


-.04 

-.OG 

.02 

-.14 

-.13 

.14 

-.02 

-.IG 

-.05 

.22 

-.03 

-.03 

.03 

-.22 

.00 

-.12 

.03 


.08 


-.01 

-.OG 

-.03 

-.10 

-.14 

.00 

.15 

.09 

.01 

.42 

-.03 


-.09 


.02 

1.00 

.04 

-.14 

.02 

.02 

-.02 

.14 


1.00 

-.01 

-.41 

-.00 

.41 

.01 



1.00 

.03 

-.00 

-.03 





1.00 

.07 


-.00 





1.00 

-.07 

.04 


1.00 .00 
1.00 


Tables 4 and 5 depict the initial and adjusted errors for the range sum rate 
biases for six stations and the survey errors for five. No survey errors were 
computed for the sixth station due to program dimension limitations. The 
established measurement accuracy goal of .003 cm/sec is met by half of the 
stations after only three days of tracking and exceeded by all thereafter. Table 5 
reveals that data from the GRAVSAT/GEOPAUSE combination can be used to 
recover survey errors with centimeter accuracy. Recovered RSS position un- 
certainties for all five stations decrease from approximately GO cm to 30 and 
20 cm alter 3, G, and 10 days of tracking, respectively. The present survey 
accuracy levels rre in the low meter range. 'I'hus significant improvements from 
the present accuracy icvlI appear possible. 

An independent analysis was conducted to d( termine the perturbations 
induced by geopotential uncertainties in the range sum rate data between a 
geosjmehronous and GRAVSAT satellite. Perturbations in the data caused by 
the differenci ' between nominal and actual geopotential values averaged more 
thim 1 mm/sec. 'I'o satisfy GRAV’SAT SEC requirements the ratio of the above 
perturbation to the perturbation caused by the unmodeled residual force should 
l)e approximately 10. Consequently, this force should cause a maximum data 
|)erturbation of 0.1 mm ec. This corresponds to at most 10 g* uneerUiintv 
for the unmodeled force. Table G reveals that alter three days of tracking an 
uncertainty of 10”*'^ g is reached and it decreases an order of magnitude for each 
additional tracking interval. 'I’hus the requirement is achieved. 

g - 980 cm sec*, the earth's equatorial otrelerotion of gravity. 
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Table 4 

Uncertainties of Range Sum Rate Measurement Biases 
After 3, 6, and 10 Days of Tracking 


Station 

A Priori 
Uncertainty 
(mm/s) 

A Posteriori 
Uncertainty (mm/s) 

3 Days 

6 Days 

10 Days 

Guam 

1.00 

.04 

.02 

.01 

Madrid 

1.00 

.04 

.02 

.01 

Rosman 

1.00 

.03 

.02 

.01 

Canberra 

1.00 

.04 

.02 

.01 

Johannesburg 

1.00 

.03 

.02 

.01 

Santiago 

1.00 

.03 

.02 

.01 


Table 5 

Station Survey Uncertainties After 3, 6, and 10 Days of Tracking 


Station 

Sl- ;y 
Component 

A Priori 
Uncertainty 

A Posteriori 
Uncertainty (cm) 


(m) 

3 Days 

6 Days 

10 Days 


E* 

10.0 

64.2 

29.1 

19.3 

Guam 

N 

10.0 

30.4 

13.1 

9.4 

V 

10.0 

21.1 

12.2 

9.5 


RSS 

17.3 

74.1 

34.2 

23.5 


E 

10.0 

28.0 

17.4 

13.7 

Madrid 

N 

10.0 

45.8 

16.9 

12.7 

14.2 

9.4 

V 

10.0 

27.7 


RSS 

17.3 

60.4 

28.1 

20.9 


E 

10.0 

32.7 

17.7 

13.8 


N 

10.0 

18.2 

9.9 

7.8 

Rosman 

16.5 

11.6 

V 

10.0 

61.3 


RSS 

17.3 

71.8 

26.1 

19.6 


E 

10.0 

50.2 

24.6 ^ 

15.7 


N 

10.0 

22.8 

16.2 

10.9 

Canberra 

V 

10.0 

26.0 

14.7 

10.2 


RSS 

17.3 

61.0 

32.9 

21.7 


E 

10.0 

34.3 

20.6 

15.6 


N 

10.0 

18.1 

10.5 

7.5 

Johannesburg 

V 

10.0 

32.4 

20.7 

13.2 


RSS 

17.3 

50.5 

31.0 

21.8 


*E, N, V are the station centered orthogonal ea«t . no* <’''.d vertica' components, respectively, ot the 
station survey error. 
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Table 6 

Uncertainty of Residual Force Upon GRAVSAT After 3, 6 and 10 

Days of Tracking 


Parameter 

A Priori 
Uncertainty 

A Posteriori 
Uncertainty 

3 Days 

6 Days 

10 Days 

Unmodeled 
Residual Force 
Upon GRAVSAT (g) 

lO"^ 

10' 

io'“ 

10"‘^ 


In (7) the coefficients of a harmonic expansion of the Goddard Earth 
Model-5 (GEM-5) geopotential field wore calibrated against actual data 
of 15® X 15° gravity anomalies. The nominal covariance matrix was 
then scaled to be consistent with the residuals. The resulting normalized un- 
certainties (a^^) from this matrix are displayed in Fig. 2 as percentages of 
Kaula's rule s^ = 10“^ /n^ , where n is the degree of the coefficient. This is an 
empirical formula which estimates the magnitude of normalized harmonic coef- 
ficients. These uncertainties increase exponentially from 5 to 60% through 
degree 12, thereafter all coefficients are uncertain to between 60% and 100% 
of their actual values. These large uncertainties are indicative of the severe 
aliasing and poor statistical independence inherent in the data from which these 
geopotential coefficients were estimated. These conditions are caused primarily 
by incomplete sampling of the gravity field. The orbit of GRAVSA'P is specifically 
designed to obtain a total sample of the field. Such a sample should allow the 
coefficients to achieve orthogonality, thereby minimizing aliasing and dependence 
between coefficients. 

Tables 7 and 8 give the improvement factors for the uncertainties of the 
sine and cosine terms of the geopotential coefficients. Improvement factors are 
the ratio of the a priori uncertainty as shown in Figure 2 divided by the uncertainty 
obtained during the simulation. Thus a factor of 100 represents an improvement 
in accuracy two orders of magnitude superior to present accuracy levels. It is 
clear there is at least one and sometimes a two order of magnitude improvement 
in the accuracy for all terms through (12, 12). Furthermori', there is only minimal 
improvement in the accuracy for most terms between different tracking intervals. 
Consequently, if accuracy improvement of geopotential terms were to become the 
objective of the GRAVSAT/GEOPAUSE mission, three days of eontinuous traeking 
apparently suffices. 
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Figure 2. Present Uncertainty of Low Order Geopotential Coefficients 
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Table 7 

Improvement Factor of the Uncertainties in Sine Term of Geopotential Coefficients After 

3, 6, and 10 Days of Tracking 
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, . A Prion Lncertainty 

Improvement 1 actor = — — — - — ^ : - 
A Posteriori Lncertainty 




















Table 8 

Improvement Factor of the Uncertainties in Cosine Term of Geopotential Coefficients After 

3, 6, and 10 Days of Tracking 
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GEOPOTENTTAL AUASING 


The geopotential field is representable by an infinite series of spherical 
harmonic terms. Estimating the field consists of solving for both the sine and 
cosine coefficients of these spherical harmonics. Since it is obviously imjwssible 
to estimate an infinite number of coefficients, this series is usually truncated 
at some term and all terms thereafter are set equal to zero or assumed perfectly 
known at some nominal value. Since neither is strictly true, the truncation 
invariably introduces an error upon the terms to be estimated which precede the 
point of truncation. This mismodeling error is called "aliasing" and is indetni 
the dominant problem in geopotential estimation. As an example of Uiis aliasing 
phenomenon, consider the following simple example. A natural phenomenon is 
modeled correctly by the (|uadratic y x^ ♦ x ♦ 1. .\ researcher, however, 
assumes that the linear equation y - ilx + b correctly mo<iels it. This assumed 
model neglects the second degree term of the correct model, in effect equating it 
to zero. Next he performs a standard least stjuares fit to the three data points 
y(0) 1, y(5) HI, and y(10) 111 using the linear model, lie determines y 1 lx - 
7.H3 Ix'st fits the data. 'Phe estimated coefficients of fit are a 11 ;mti b -7.HH 
whereas the correct values are a b 1. 'Ihus neglecting the second degree 
term in tlie correct motiel has seriously degraiieil the (]uality of the estimates 
from tlie least scpiares fit. This degradation is an example of the "aliasing". 

This same aliasing exists on a much larger scale in geo|H)tential coefficient 
estimation. The neglect of the uncertainties of the high degree and order unesti- 
mated coefficients aliases the error of the remaining estimated coefficients. A 
mathematically rigorous description of this aliasing and orthogonality properties 
of geopotiuitial coefficients is ilevelopi'd by Argentiero and Garza-Uobli's (Hef. S). 


IXiring a eo\ariaiuu* analysis each parametiu' uiuler investigation is placed 
in 1‘itlier an "cstimati'd" or "i-onsid«.*r" catt'gory. 'I’lu* fornuu* contains tliosc 
paramctiu’s whieli ordinarily wouhl solvi'd for during an actual orl)it determi- 
nation process, c.g. satellite state or instrument l)iasi*s. Tlu' latter includes 
parameters whieli winilil remain unestimatial tluring such a process Inil wliusi- 
uneertaintit's are takiui into account in computing the covariance matrix of 
the process. Measurenuuit timing ami atmosplierii- i i'li'aetion ari* typiiuil 
eonsidi'r paranu'ters. The software (Ib'f. used »luring this study I'on- 

veniently separates the variance ol each ''estimated" parameter into a sei ies 
of variaiu-es due to tlu' individual "eonsidei" pai ameters and tlu* vai ianei' due 
to the data noisi*. The i*rror i*ontril)ution ol a spi'iulii’ "consider" parameter 
to till' error of an "estimated" parameter is called its "aliasing." 


/ 


t 


;\s mentioned previously, results to ilate neglect the aliasing phenomenon 
introduced by the uiu-ertainties of the unadjusteil parameters and thus lend to be 
somewhat optimistic. To oblr.in some quantitative measure ol this aliasing, the 


covariance analysis was reconsidered using two different estimation strategies. 

The first strategy attempts to determine the aliasing impact upon the uncertainty 
of a centrally located adjusted geopotential term of low degree and order due to the 
uncertainties of the surrounding unadjusted terms. Consequently, the uncertainties 
of both satellite epoch states, measurement biases, station survey components, un- 
modeled residual force, and an arbitrarily selected geopotential term, C(5, 2), were 
placed in the "estimated" category. The uncertainties of the remaining 164 
geopotential terms surrounding C(5, 2) were placed in the "consider" mode. Such 
a strategy reveals, at least for one particular term, the degrading influence of 
the remaining terms and would identify those particular ones with the greatest 
aliasing. Alias maps for all terms should be estimated for the total aliasing 
effect. However, an alias map for one term gives at least some insight to the 
aliasing problem. The aliasing effect of the uncertainties of the unadjusted terms 
upon the error of C(5,2) is displayed in Table 9 as an alias map for C(5,2). The 
unnormalized a priori uncertainty (aj,) of C(5,2) as determined from Figure 2 is 
10‘*. The adjusted uncertainty (5^) was 10“’. In general, two results are immediately 
obvious. First, the uncertainties of all the surrounding terms alias the adjusted 
error of C(5, 2) to varying degrees. Second there is wide variation in the magnitude 
of this aliasing. This second factor is of special importance since it indicates 
which particular geopotential terms should be estimated concurrently with C(5, 2) 
to minimize the overall degrading influence of aliasing. At least to the degree and 
order considered in the alias map, an interesting pattern immediately emerges. 

The surrounding terms of odd degree and even order to order 4 exhibit markedly 
higher aliasing upon C(5, 2) than the remaining terms. The contributions of these 
particular terms are boxed twice within Table 9. The magnitudes of the aliasing 
for these terms is greatest for the lowest order and decreases rapidly within a 
degree reaching approximately 2 or 3 x 10“‘° HSS^*^ aliasing at order 4. Within 
an order, the aliasing increases slightly with increasing degree. Indeed, term 
(11, 0) alone contributes 182.4 x 10”‘° in an HSS sense to the adjusted uncertainty 
of 10"’ for C(5, 2). Thus the odd degree zonals are seen to contribute the maximum 
aliasing upon C(5, 2) and therefore should be estimated simultaneously with C(5, 2). 
All the remaining terms contribute much less aliasing and could remain in the 
consider mode. 

The second strategy investigates the extent and magnitude of aliasing induced 
by an unestimated high degrt*e and order geopotential term uix)n terms of lower 
degree and order. The covariance analysis was reconsidered with only the un- 
certainty of C(12, 0) in the consider mode :md all remaining uncertainties 
estimated. C(12, (») has an initial uncertainty of 10"'^ . The information obtained 
by such a strategy is important because it determines how much of tne 


(1) RSS of sine ond cosine components of the geopotenfiol term. 


Hi 



Table 9 

Alias Map for C(5, 2) 

RSS Aliasing Contributions x io‘® to the Uncertainty of 
C(5, 2) Due to the Uncertainties of Surrounding Unestimated 
Coefficients After 10 Days of Tracking 

Order (M) 


0 1 2 34567 89 10 11 12 


2 

3 

4 


5 


Z 6 


(U 

^ 7 

hfl ' 


a 


8 

9 

10 
11 
12 



uncertainty of an adjusted term is in fact due to the aliasing by the uncertainty 
of an arbitrarily selected unadjusted term. 

The results of this strategy are displayed in lx)th Figure 3 and Table 10 
as percentage aliasing of the estimated uncertainty. One would intuitively ex|)ect 
that with such a dense and global distribution of data as obtained by (iHAVSAT/ 
GEOPAUSE, the aliasing would be negligible. It is seen, however, that even 
with such a good distribution of data the aliasing problem is quite severe. 
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PERCENTAGE 

aliasing 



Figure 3. Histogram of Percentage Aliasing in Uncertainty of Estimated 
Cosine Coefficients Due To Unestimated Uncertainty of C(12, 6) 



Table 10 

Percentage Aliasing in Uncertainty of Estimated Ck^sine 
Coefficients Due to Unestimatcd Uncertainty of C(12,6) 


Order (M) 



0 

1 

2 

3 

4 

5 

G 

7 

8 

9 

10 

11 

12 

2 

.01 

.53 

.14 










3 

.08 

D 

.8G 

.20 









4 

.02 

2.1 

.47 

.43 

1 

1.6 








5 

.G2 

1.9 

.08 

.2G 

.24 

.63 







6 

.13 

2.3 

.10 

.73 

2.1 

.01 

13.4 






7 

.G4 

3.0 


.71 

1.8 


2.0 

.17 





8 

.72 

.29 

.27 

.59 

1.1 

3.8 

19.9 

.79 

1.2 




9 

.08 


.25 

.18 

1.1 

5.3 

1.1 

.55 

1.8 

.40 



10 

.G3 

? ■ — 
.07 

.45 

1.2 

7.0 

1 .3 

71.8 

.94 

.98 

1.2 

.14 


11 

.10 

3.8 

1.4 

.93 

.25 

8.1 

.26 

.71 

1.2 

..34 

.3.3 

.17 


12 

.34 

.13 

.IG 

.57 

2.3 

1.7 

0 

lO"*-* 

.95 

2.3 

2.G 

.12 

.13 

.11 


The unadjusted error ofC(12,6)is seen to alias the estimated uncertainty of all 
terms to some degree. 'ITie uncertainties of seven of these terms: C(G, G), 
C(8, G), C(9, 1), C(9,5), C(10, 4), C(10, G), and C(ll, experience the greatest 
impact with percentage aliasing ranging from 5 thmugh 72'i. 'Hie uncertainty 
for C(10, 6) stands out with the highest percentage aliasing of 72['i . ITie uncer- 
tainties of all the remaining terms experience an aliasing effect lx?low 3'". 
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CONCLUSIONS 


This report has described a covariance analysis of the GRAVSAT/GEOPAUSE 
mission in which the uncertainties of approximately 200 parameters were simul- 
taneously estimated from the range sum rate data between these two satellites. 

The adjusted GRAVSAT epoch state errors are seen to be sensitive to the amount 
of data w’hereas the GEOPAUSE errors remain relatively insensitive. The actual 
magnitude of these errors indicate that data from the high GEOPAUSE spacecraft 
to a lower satellite can successfully be processed to estimate both orbits. In 
addition RSS position errors ranging from GO to 20 cm, depending upon the amount 
of data, can also be extracted from such data. These arcuracy levels represent 
improvements over the present 8.5m RSS accuracy level. The uncertainties of 
the geopotential field to (12, 12) were recovered to an accuracy one to two orders 
of magnitude superior to present estimates regardless of the amount of data. 
Consequently, if geopotential recovery is of prime interest, 3 days of tracking 
suffices. Geopotential aliasing is the degrading influence upon the uncertainties 
of estimated parameters due to the errors in the unestimated quantities. An 
alias map for an arbitrarily selected term, C(5, 2), revealed that the uncertainties 
of odd degree zonals contribute the maximum aliasing to the uncertainty of that 
term. The percentage aliasing in the estimated uncertainties ol coefficients due 
to the error of a centrally located term of high degree and order term, C(12, G), 
was also determined. The estimated uncertainties of most terms contained 
negligible, i.e. less than 5^^, aliasing. However seven terms, ranging from degree 
G through 11, contained significantly higher aliasing. 
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